The initial singularity of ultrastiff perfect fluid spacetimes 

without symmetries 



J. Mark Heinzle^* and Patrik Sandin^^ 

^ University of Vienna, Faculty of Physics, Gravitational Physics, 
Boltzmanngasse 5, 1090 Vienna, Austria 

2 

Department of Physics, University of Karlstad, 
^Nj ' S-651 88 Karlstad, Sweden 

o 



January 20, 2013 



Abstract 

CJ 

We consider the Einstein equations coupled to an ultrastiiT perfect fluid and prove the existence 
of a family of solutions with an initial singularity whose structure is that of explicit isotropic models. 
f^£y This family of solutions is 'generic' in the sense that it depends on as many free functions as a general 

solution, i.e., without imposing any symmetry assumptions, of the Einstein-Euler equations. The 
^ , method we use is a that of a Fuchsian reduction. 

> 
m 

: 1 Introduction 

1^ In general relativity, spacetime singularities can occur under very general conditions, see [1]. However, 

\ mathematical theorems detailing the structure of these singularities usually assume a great deal of sym- 

metry; either spatial homogeneity is assumed, or as in the Gowdy and cases, an Abelian two-parameter 
isometry group acting on spatial sections. A potentially powerful method to obtain theorems concerning 
non- homogeneous models is the method of Fuchsian reduction [2], which is a rigorous attempt to 'expand' 
^ solutions in the neighborhood of a singularity. A particular way (but by no means the only one [2]) to 

r> ■ obtain the leading order behavior of solutions toward a singularity is to formulate an 'asymptotic system' 

\ of equations that is derived from the Einstein equations by dropping terms, e.g., spatial derivatives, which 

" ■ ■ are believed to be unimportant for the leading asymptotic dynamics. If this 'asymptotic system' truly 

captures the leading divergent behavior of solutions of the Einstein equations toward a singularity, then 
the remaining part of each solution obeys equations that are of so-called Fuchsian type. In this case, it is 
possible to prove the existence of a general class of solutions of the Einstein equations whose asymptotic 
behavior is governed by the behavior of solutions of the asymptotic system. This method of Fuchsian 
reduction was applied, e.g., in [3] to the initial singularity of cosmological models with a stiff fluid or mass- 
less scalar field in 3 -I- 1 dimensional spacetimes without symmetries, and to the Einstein-dilaton-p-form 
system in higher dimensions [4]. 

The structure of spacetime singularities is particularly interesting in the context of the Weyl curvature 
hypothesis, where the focus is on isotropic singularities [5]. Spacetimes with such singularities are not 
generic among solutions of the Einstein- matter equations, i.e., an arbitrary matter/curvature distribution 
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at late times does not develop an isotropic initial singularity when traced back in time. However, this 
is expected to only hold for matter sources where the ratio of pressure to energy density is less than or 
equal to one; if this ratio is greater than one, isotropic singularities are expected to be generic. In the 
case of a perfect fluid matter source, such a matter model is called an ultrastiff perfect fluid. 

The interest in ultrastiff fluids has come from attempts to obtain non-singular descriptions of spacetimes 
representing collapsing universes bouncing outward before a singularity occurs. This in the case in the 
ecpyrotic scenario or the cyclic brane world models [6] , where the ultrastiff fluid is an effective description 
of a scalar field with negative potential energy during the contraction phase. The scalar field is in turn 
a representation of the size of a compactified dimension and the potential a representation of the force 
between two different brane worlds. These models have led to an interest in ultrastiff cosmologies in 
general [7] and especially in their structure in the final stages of collapse, or equivalently, initial stages 
of expansion in an expanding universe. The structure of the initial singularity in cosmological models 
with an ultrastiff perfect fluid has been studied under the assumption that the singularity is 'local' [8, 9], 
which means that the time evolution at adjacent spatial points decouples. Under this a priori assumption 
of locality it has been demonstrated that the past asymptotic structure is described by an isotropic 
solution [8, 10] of the Friedmann- Robertson- Walker type, i.e., that anisotropics become small compared 
to the overall expansion of the universe close to the singularity. The situation is similar for cosmological 
models with several non-interacting perfect fluids, of which at least one possesses an ultrastiff equation 
of state [11]. 

In this paper we investigate the past asymptotic dynamics of cosmological spacetimes with ultrastiff 
perfect fluids. We do not impose any assumptions on the symmetry of the spacetimes and we do not 
make any a priori assumptions on the past asymptotic behavior (like asymptotic locality) . We prove the 
existence of a family of solutions that converge to isotropic models toward the initial singularity; the past 
asymptotic behavior of these solution is characterized in detail, see Theorem 4.1. This family of solutions 
we construct is 'generic' in the sense that it depends on as many free functions as a general solution of 
the initial value problem connected with the Einstein- Euler equations. 

The paper is organized as follows. In section 2 we give the Einstein-Euler equations in a Gaussian 
coordinate system; in section 3 we describe the asymptotic system and its role in the Fuchsian reduction. 
Section 3.1 is devoted to the special case of a perfect fluid with equation of state p = 3/i, for which 
the asymptotic system admits explicit solutions. We use these solutions to set up a Fuchsian system of 
equations; this is done in analogy with the analysis of the stiff fluid case [3]. In section 3.2 we give the 
treatment of the general ultrastiff fluid case, which is less explicit but follows the same principles. The 
main result is then stated in section 4, while section 5 gives a derivation of the reduced equations that 
are used to prove the main result in section 6. 

2 The Einstein-Euler equations 

Let be a four-dimensional Lorentzian manifold with metric '^gap (with a, (3 ~ 0, 1, 2, 3). We consider 
the Einstein-Euler equations for self-gravitating perfect fluids, i.e.. 



where we assume a linear equation of state p = w/i, w = const, relating the energy density fi (as 
measured in the fluid's rest frame) and the pressure p of the fluid; it" is the four-velocity of the fluid. 
We use geometrized units, i.e., c = 1 and SttG = 1, where c is the speed of light and G the gravitational 
constant. 

We consider spacetimes that are diffeomorphic to R x M, where M is three-dimensional, and metrics of 




(la) 
(lb) 



V^T"/^ = , 
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the form 

-dt" + gabL0''u\ (2) 

where gab = gab{i) denotes a one-parameter family of Riemannian metrics, which are naturally identified 
with the metrics induced on t = const hypersurfaces; the coframe {cj^,cj^,a;^} on M is arbitrary. The 
range of Latin indices a,b, . . . is 1, 2, 3. 

The Einstein-Euler equations become a first order system of equations with constraints when the extrinsic 
curvature kab — —\dtgab of the t = const hypersurfaces is used (3 + 1 split). For our purposes it is 
preferable to split kab into its trace ('mean curvature') and (the negative of) its traceless part, i.e., 

kab = -(yab + \ (tr k) gab ■ (3) 

The tensor cjab coincides with the rate of shear tensor of the (geodesic) congruence orthogonal to the 
t = const hypersurfaces. Then the evolution equations are 

dtgab = '2.gac{(T\ - i(trfc) 6\) , (4a) 

dta\ = (tr k)a\ - {R\ - ^R5\) + {S\ - i(tr S)S\) , (4b) 

9t(trfc) = i? + (trfc)2 + KtrS) - fp, (4c) 

and the constraint equations read 

R -a\a' a + lit! k)^ ^2p, (5a) 

-Vaa", -fV,(trfc)=jfe, (5b) 

where 

p ^ fl{l + {1 + ■w)u"Ua) , jb = p{l + tv){l + UaU^y^^ Ub , Sab ^ Pt{{i + w)UaUb + Wgab) ■ (6) 

The quantity p is the energy density, ja the current measured by the observer associated with the geodesic 
congruence. Equations (4) and (5), together with (6), correspond to (la). The Euler equations (lb) read 

dtp, — (1 + w){tTk)p = —(1 + w) v?dtp + 2pu'^dtUa — pcjabu'^u'' — |/i(trfc)M^ 



+ Vl + u"VaA* + y"(l + u^y^/^Uau''VbU°' + pVl + VaU" , (7a) 

dtUa + W (tr k) Ua + Y^/^^^VajU = -p^^Ua [dtP " (1 + w) p{iv k)] -Ua{l + U^)~^ [u^ OfUb - (Jbcu''u''] 

+ \ua{l + U'y\tvk)u' - [(1 + u2)-l/2 _ 1]_^^-1V,^ 
- (1 + u2)-l/2 [^-^UaU^VbP + UaVbu'' + V^V bUa] , (7b) 

where we use the abbreviation = UaU'^- 

3 The asymptotic system 

Alongside the Einstein-Euler equations we consider another system of equations for the same set of 
variables, which we denote by ( °gab,tT °cr"jj, %a) in the present context. This asymptotic system 
consists of evolution equations for ( "^afaitr "fc, 'V"^), 

dt %b = 2 °gac ( "a^ - i (tr °k)6%) , (8a) 
dt °<J\ = (tr Ofc) , (8b) 
dt (tr O/fc) = (tr 0fc)2 + |(u; - 1) 0^ , (8c) 
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and evolution equations for the matter variables ( ^^ta), 

c)tV-(l+w)(tr °/fc)V = 0, (9a) 

dt V + w (tr Ofc) °Ua + V"' °Va V = . (9b) 
In addition, we consider the constraint equations 

2V, (10a) 

% , (10b) 

where °p = and ^ja = (1 + w) %a. The symbol °Va denotes the covariant derivative associated 
with the metric '^gab- 
Note that the evolution equations (8) and (9) preserve the constraints (10). A computation shows that 

dt{~ V, V, + |(tr \f - 2 °p) = 2(tr %)(- V, V, + |(tr "kf -2"p), 

dt{- V. V, - |V,(tr Ofc) - %) = (tr "fc)(- V, - I V,(tr %) - \) 

-i'V,(-V,V, + |(tr°fc)2-2°p), 

hence if the constraints are satisfied initially, they are satisfied for all times. 

The asymptotic system (8)-(10) is a system of partial differential equations, where, however, the spatial 
derivatives merely enter in a passive manner, through the decoupled equations (9b) and (10b); hence 
in spirit, (8)-(10) are ODEs. It is clear that the asymptotic system is obtained from the Einstein- 
Euler system by dropping (a large number of) terms. The derivation of the asymptotic system from the 
Einstein-Euler system is intimately connected with considerations involving the concepts of asymptotic 
silence and locality; we refer to [9] and references therein. The aim of these considerations is to find 
a simple system of equations that reproduces the dynamics of the full Einstein-matter equations in the 
asymptotic limit toward a 'generic spacelike singularity'. In other words, a simple ('asymptotic') system 
is sought that governs the asymptotic behavior toward a spacelike singularity of (typical) solutions of 
the Einstein-matter equation in the sense that each (typical) solution of the Einstein-matter equation is 
approximated, to a certain order, by a solution of the asymptotic system. 

Whether the asymptotic system given in (8)-(10) indeed governs the asymptotic behavior of solutions 
of the full Einstein-Euler equations, is a non-trivial question. The considerations of [9] strongly indicate 
that a necessary condition is that the equation of state parameter satisfies w > I. (In the case w < 1, the 
conjectured asymptotic system is considerably more complicated, because it must capture the conjectured 
oscillatory behavior, 'Mixmaster behavior', of solutions toward the singularity.) The case of a stiff fluid, 
i.e. w = I, has been treated in [3]: The asymptotic system given in [3] is shown to govern the asymptotic 
behavior of solutions of the Einstein-Euler equations with a stiff fluid. (Note, however, the comment in 
the symmetry subsection of section 3.2.) 

In this paper we consider ultrastiff fluids, which are characterized by w > 1. We prove that the asymptotic 
behavior of solutions (more specifically, of a full set of solutions in the sense of function counting) of the 
Einstein-Euler equations with an ultrastiff stiff fiuid is governed by the asymptotic system (8)-(10); in 
addition, we give detailed information about the order of the approximation. 

3.1 A particular (ly simple) case of ultrastiff fluids: w = 3 

We begin by solving the asymptotic system with w ~ 3. The reason for this choice is one of convenience, 
since in the special case w = 3, the solutions of the asymptotic system take a simple and explicit form 
(which enables us to stress the analogies of the present analysis with the analysis of [3]). The general 
(but less explicit) case w > 1 is treated in subsection 3.2. 



-VfcV, + |(tr°fc)2 = 
-VaV,-|V,(tr%)== 
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We are interested in solutions with '^^ > and tr \ < 0, i.e., 'expanding' models with positive energy 
density. (Note that "/i remains positive if it is positive initially, see (9a); likewise, if tr °fc is negative at 
t = to, it remains negative for all t < to because 9t(tr °fc) > 0, see (8c).) The basic observation is that 

dt{\/i^±tr °k) ^±{y^±tr °kf. (11) 

The solutions of (11) with the minus sign are 

v/3^-tr"fc=— , (12) 
t + tplx) 

with t > ~Tp(x), where ip = ip{x) is an arbitrary function of the spatial coordinates which we collectively 
denote by x. (The trivial solution is excluded by our assumptions.) The 'initial singularity' is signalled 
by the divergence of (12) at t We wish this initial singularity to occur at i = 0. To simplify 

further computations we thus choose our initial hypersurface to be a surface where ip = const. This is 
not a restriction of the generality of the result since it is always possible to find such a surface between 
any given initial hypersurface and the singularity. The constant may be set to zero by a shift in the time 
coordinate. 

The solutions of (11) with the plus sign are 

^/3^ + tr "fc = 5— and + tr °fc = , (13) 

t + (j){x) 

where (f> — (t>{x) > is an arbitrary function of the spatial coordinates. (The function is positive, since 
2-\/3 = t^^ ^ {t + 4')^^ > 0.) The trivial solution in (13) can be thought of as arising as a special case 
by taking the limit — > oo. We conclude that 

Integration of the remaining evolution equations of (8) and (9) then yields 

02.a 

Vf, = , (14b) 

V= '{la{t{l+r't)f^-5(b-^{ya(b)t^l + lr't), (14c) 



"gac{t (1 + 0-ii))'/'(v/F^ + v/1 + (14d) 



The equations for the metric and its inverse use the matrix exponential z-^ := e^''^s^)A^ which is well 
defined for any positive scalar z and square matrix A. Note that {^/ (f>~^t + ^J'i+ljpH)^ in the 

limit — >■ oo, hence ^gab = °5act^^'^ cxp(4v^ "(t^^) in that special case. 

The contraints (10) give restrictions on the free functions of the spatial variables: </>, ^o""^, %ah^ °tta- We 
have "(T^f, ^o-^^ = 2(/)/3 from (10a). The constraint (10b) can be used to express, e.g., in terms of "o-'^^, 
"^^afa, and 0. For the Fuchsian analysis of Section 6 we will assume that the free functions are analytic in 
the spatial variables x. 

3.2 Ultrastiff fluids with general w '> \ 

In the case of a general ultrastiff fluid with an arbitrary value of w > 1 the solutions of the asymptotic 
system, which are given by (14) in the case w = 3, cannot be given in explicit form. To derive the 
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solutions of the asymptotic system with w 7^ 3 we resort to power series; the ideal framework is the 
Hubble-normalized dynamical systems approach. 

Consider the asymptotic system (8)-(9). The equations (8a) for ^gab and (9b) for ^Ua, together with (10b), 
decouple from this system, i.e., the asymptotic system reduces to equations for °cr%, tr and which 
are subject to the constraint (10a). We write these equations in the standard Hubble-normalized variables, 
i.e., we introduce 

and an adapted time variable d^- = H^^dt- In these variables we obtain a decoupled evolution equation 
for H, 



H=-\ tr \ and T."^ = , Q = — ^ , (15) 



drH = -^H{w + 1-{w-1)J:^) , (16a) 

and an evolution equation for E^j^, 

dr^\= l{w-l)E\{l~E'), (16b) 

where E^ — iE°^S^^. The constraint simplifies to = 1 — S'^ (which makes the evolution equation for 
51 redundant). 

Equation (16b) can be solved straightforwardly; we obtain 

^ " - ^ " e2(--i)- , (17a) 



where E"^ are functions of the spatial variables. Furthermore, (16) yields 



H^l - E2 = nyi - E2 e-^^'^'+i)^ , (17b) 

where H depends on the spatial variables. Consequently, 

°(T°f, = Ht^'f^e'^^ , (18a) 



tr °fc = ~iH = -3H y^l-E2 e-i(i+"')" (l + e^^'""!)") , (18b) 

and "/i = 3i?2(l — g-3(!ii+i)T^ Accordingly, we obtain explicit expression for the solutions of the 
asymptotic system in the time variable r. 

In order to express (18) in terms of cosmological time t we must integrate the relation dt = H~^dT where 
H is given by (18b). For large negative r we can expand the square root in a power series in ^^i^-^)'^ ^ 
and integrate it term by term. Requiring synchronous cosmological time, where t — Q corresponds to the 
initial singularity, we obtain 



2e2 

t = 



^A{l + w)H^Jl-t^ 

\ 2{3w - 1) 1 - E2 
as T — ^ — 00. Inverting this relation results in 



3(1 +W) / T."^ \ 2 6(«,-l)r , (J(„9{w-1)T\ 



g2(l+»)r ^ 3 (1 ^ nj)HJl - E2 i X 



/ l+W/3(l + w)\2(."-l)/(l+«') J^2jj2{w-l)/{l + w) 2(w-l) 4(^-1) 

X l + ^^(-^ -] , t +0(t 1+"' 

I 2 \ 2 J (3w- 1)(1 -E2)2/(i+-) ^ 
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in the limit t — > 0. A calculation then shows that 

= . / . \ -2/(l+Mi) „ , 2 , 2(t»-l) , 

2 . 2(w-\) , 

EE V^i i+t" (l + 0(< 1+"- )) , (19a) 



tr \ = -3H = 1 



1 /3 1 + S2if2(^-i)/{i+-) 2(^ 4(^-1) 

- '-] ^ , 1 1+"' +0(t i+t" ) 

1 V 2 / (1 - s2)2/(i+.u>) V ; 



1 + w \ 3ti!-lV 2 / (1-^2) 

2 / 2(m;~1) 4(-m-1) \ 

ee; 1 + 0-4 1+"- +0(t 1+"' ) , (19b) 

1 + w \ ' ) 

where we have introduced the spatial functions "^(T^^ and In addition, (8a) and (9) yield 

4 / 7/1+1 w-\ 2(m-l) 



(19c) 



4 1 /?/; + ! 2(«.-l) 4(u;-l) \ 

"^^ = 17?^*"' i-!li±l0-itT+^+O(iTT:r) , (i9d) 
3(1 + wy \ w — \ J 

2w w(w + 1) 2m tti-1 2w 2(id-1) 

V = °Uat^+^ - -^<i)^'^\/a(t) t^+^t^+w +t^+wO(t l+J" ) , (19e) 

(w — \y ^ ' 

and requiring that ^g°^ is the inverse of ^gah gives 

2(m)-1) 



4 / 7/1+1 w~\ 2(m)-1) \ 

V'= 0^»^t3(l + u;) (^b^_2^p^ +0(<^+^) j 



(19f) 



The expressions (19) replace (14) in the context of the case of general w > 1. It is not difficult to convince 
oneself that (19) reduces to (14) in the special case w — i. 



4 Main result 

Stated in an informal way, the main result we will prove is that ('typical') solutions of the Einstein- 
Euler equations 'look like' the solutions (19) of the asymptotic system in the asymptotic limit toward the 
singularity. The rigorous statement is the following. 

Theorem 4.1. Let = {^gabM °k, ^(t'^jj, ^fi, ^Ua) be a solution of the asymptotic system as given 
by (19). Then there exists a unique solution W = (^gab,^^ k,a°'ij, iJ,,Ua) of the Einstein- Euler equations, 
whose asymptotic behavior, to leading order, is that of "W. Specifically, there exists a set of quantities 
U = (7°t, "J^b, V, Va) that vanish in the limit i — > 0, together with their spatial derivatives, and an 
arbitrarily small e > 0, such that 

0„ , j.2-4/3(H--u))-2e 0„ (on^\ 
gab ^ gab + t ' ' gad b^ i^Oa) 

tr = tr Ofc + <2-io/3(i+«,)-e ^ ^ (20b) 

= + il-4/3(l+-)-e/2 ^a^ ^ (20c) 

/i= V (20d) 

Ua^ V+i^-^°/'(^ + -^-'^«a. (20e) 

The inverse metric is given by g"'' ~ ^g"'' + i2-4/3(i+iu)-2e ^g'^'^ , where 7°^ shares the properties ofU. 
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Remark. The number of free functions (of the spatial variables) in the set of solutions of the asymptotic 
system is the same as the number of free functions that determine initial data for the Einstein-Euler 
equations. Therefore, in the sense of function counting, Theorem 4.1 shows that 'typical' ('generic') 
solutions of the Einstein-Euler system behave according to (20). 



5 The reduced system 

The proof of Theorem 4.1 is based on techniques in connection with Fuchsian equations [2, 3, 12]. Our 
aim is to construct, from the Einstein-Euler evolution equations (4) and (7) and the asymptotic sys- 
tem (8)-(10), a system of (generalized) Fuchsian type, i.e. a system of the form 

tdtU + AU = Tit,x,U,dM) + g{t,x,U)tdtU (21) 

where U denotes the collection of the variables of Theorem 4.1 (and additional variables). In this context, 
A is required to be a time-independent matrix that depends analytically on the spatial coordinates, and J- 
and g are functions that are continuous in t > and analytic in the other variables. SufBcient conditions 
that the system (21) be a (generalized) Fuchsian system are that (i) the matrix A is the direct sum of a 
submatrix whose eigenvalues have positive real parts and the zero matrix, and that (ii) the functions J- 
and g vanish like a positive power of t as t — > 0. If this is the case, then there exists a unique solution 
U of (21) such that U and its spatial derivatives vanish with t; see [2, 3]. Our aim is to construct a 
system (21) and to show that the 'Fuchsian conditions' are met. 

Let "VV = ( ^gab, °5°'',tr °cr°j, ^Ua) be a solution of the asymptotic system as given by (19). Let 
us (re)define U asU — (7";,, 7"j,, k, Va, ^°'bc^ ^"hc)- We make the ansatz 

5afc= Vb + i"^ VcT'b, (22a) 

gab ^ O^ab ^ -a^ OgCb ^ (22b) 

tr/c = tr °fc + t"''+^"WT) (22c) 

CT^fc = "(T°b + r^'+^-H^ , (22d) 

M = V + t""^^^(^ V , (22e) 

Ua = V + i2-/(i+-)+". „^ , (22f) 

7V = ^~"^A^,, (22g) 

7V = t-^^A?,, (22h) 

where (a-y, a-y, a?, Q^i/, an, Q;a, Q^a) are constants. 

The procedure is (almost) straightforward: Imposing the Einstein-Euler equations yields a system of 
equations for the variables U. The plan is then to find a range of the free constants such that this system 
is Fuchsian, which in turn implies the existence of a unique solutionis with the boundary condition — >■ 
(and dJA -> 0) as t 0. 

However, a slight modification of the Einstein-Euler equations is needed for this purpose; cf. [3]. The 
problem concerns the symmetry of the metric tensor gab- Obviously, for an initial value problem, the issue 
of symmetry does not arise: The Einstein evolution equations propagate the symmetry of the metric (and 
the second fundamental form). In the present context, however, the solution U of the Fuchsian system 
we intend to obtain, is not determined by its initial value but by requiring the boundary condition — > 
(i — >■ 0). It is thus not clear, a priori, whether the tensor g^b of (22a) will be a symmetric tensor, since 
we do not have control over the properties of 7°^ . This is turn suggests that special care is needed in 
connection with the curvature terms in the Einstein equations: In (4) and (5) we replace the curvature 
terms TC^ ^"^^ ^ by the corresponding terms that are defind from the symmetric part of the metric 
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^dab = 9(ahy- i-G-, R°'i, ^^"f, = R°'b[^9cd\- Thereby we modify the Emstein-Euler equations in such 
a way that the equations make sense for general non-symmetric tensors but still reduce to the original 
equations in case of symmetric data. In addition, we must establish conventions concerning raising and 
lowering of indices with a metric tensor that need not be symmetric. We define g""^ as the unique inverse 
oi gab, gacg'^^ = and use the convention that indices on tensors are raised and lowered with the second 
index of g°'^ and gab, respectively. As a consequence, index raising and lowering are inverse operations, 
as is the case for symmetric metrics. 

Having prepared the Einstein-Euler equations in this way we can insert the ansatz (22) into (4) and (7) 
without difficulty, and we obtain the reduced system, a system of equations for the variables lA. (By 
construction, this system is well-defined irrespective of the symmetry properties of the variables lA. We 
will return to the question of symmetry of the unique solution we subsequently construct in a separate 
subsection.) 

The reduced system reads 



(23a) 
(23b) 



X + e-' <;\rb - i t"''^"^- « {s\ + e-' rb)), 

tdtK + (1 + a^)K - i(l + Sw;)^"--"'' ly = 2^-"^-"^ ^\ V, + [l + ^ tr %] k 

+ (l + u;)t''/3(i+"')-"» °^M2 + (l + it;)r''-"«u2^, (23c) 

tdt^\ + (1 + - Oi,Y\ = 0(7% K + i2(l+3t«)/3(l+«,)+a.^a^ ^ 

+ _|_ [1 + 1±H t tr Ofc] ^\ - t4/3(l+-)-ha. (S^a^ _ fRb\) 
+ (1 + ti;)t4/3(l+-)+a, O^^^a^^ „ 1 ^2^a^) 

+ (1 + w)t"''+"^- (w^tb - \ u^i^by, (23d) 

tdtV + (2 - + 0.,)V - t"»-"^ K ^ 2 [1 + i±H t tr "fc] + (1 + ^)i2(l+3»)/3(l+»)+a. 

-3(TW^"''^"il-t(l + -)'"^^']'^ 

+ ^1+4/3(1+'^)-". [r.h.s. of (7a)] , (23e) 

tdtVa + a^Va = - ^[1 + i±Ht tr \\ Va ~ W + K ^a " wt^^^+^-^/^^^+-^+"^ KVa 

- _J^( 0u)-lt-l+2/3(l+-)+"^— "V„i. + j^i )-li-H-2/3(l+«>)+a.-a. x 

X ( V + i-*/3(l + "'+"-Z^)"'(Va V + VaJ.) J. 

-f il-2u,/(l+«,)-a„ [j. 1^ ^ (23f) 

t9tAg, - t°^t9t7%,c + aAt"^7%,c , (23g) 
tdt'X^ = eHdtTb,c + aAt"^7%,c. (23h) 

For brevity we do not write out explicitly the terms on the right hand sides in equations (23e) and (23f). 
The expressions in question are straightforwardly obtained by inserting the ansatz (22) into the r.h.s. of (7). 
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6 Fuchsian analysis 



6.1 The reduced system as a Fuchsian system 

In this section we perform the proof of Theorem 4.1. The first part of the proof of is to show that there 
exist appropriate choices of the constants (a-y, a-y, a^, a?, cn^, a^, ax, ax) such that the reduced system (23) 
becomes a Fuchsian system. 

The first step is to consider the Lh.s. of (23). As a prerequisite, the Lh.s. must be of the form tdtU + AU, 
where W denotes the collection of variables, i.e., W = (7°^, j^, k, (;\, Ua,A°^g, A°{,^). This requires 
setting 

= a^ =■■ P ■ 

To establish the conditions on the coefficient matrix A that are required for the system to be Fuchsian, 
cf. [3, Section 4], we show that A is the direct sum of a positive definite matrix and the zero matrix, 
provided that the constants {a^, a-y, a^i^ ct<;i C(i^[~ I3],ay,ax, ax) are chosen appropriately. 

Obviously, the only non-diagonal part of A, associated with the variables {k,i'), corresponds to the 
submatrix 

1 + 13 -i(l + 3w) 

-|(l + u.)-i 2- l{l+w)-^ + P 

The eigenvalues of this matrix are /3 and 3 — ^{1 + w)^^ + f3, which are positive if and only if /3 > . There- 
fore, the conditions ensuring positive definiteness of the submatrix associated with (7"f,i7"6i '^,Va) 
are 

2 

a) > b) > c) a,^ = a^ ^ > d) < 1 4- — e) a„ > . (24) 

0(1 + w) 

In the second step we turn to the r.h.s. of (23), which is of the form T + Q tdtU, where T depends on t, 
the spatial variables, and U and its first spatial derivatives, and Q = Q{t,x,L{). We need to show that 
J" and g e 0(t*) for some positive number 5, if the constants (a^,a-y,a„[= /3],a^,a^[= 13], ay, ax, ax) 
are chosen appropriately. We use the notation f{t,x) £ 0(^g{t)) to compare the asymptotic behavior of 
a function f{t,x) and a function g{t) near t = 0, if f{t,x) = 0(g(t)) as t — > uniformly on compacts 
subsets of the spatial variables x. 

Equations (23a) and (23b) imply the conditions 

2(14- 3w) _ 2(1 3m;) 

a) "'^+"^<^(TT^ b) a, + a,<^^^. (25) 

To treat (23c) we use that differs from ^u^ by a function of time that goes to zero as t — > faster 
than itself, cf. (22b), (22f), and (24b, e). The leading order is thus obtained by replacing by 
in (23c). We get the conditions 



a 



) ac+/3<l-TT- 2ac+/?<2- ^ c) /? < 2 - . (26) 

1 + w 3(1-1- uO 1 + w 



Analogously, (23d) leads to 



/ 2x1 — w , / 4x 1 — w , ^ 

a a^+p >-(!-- ) = - b ac>-(2-- )=2- ; 27 

' ^ ^ l + w' 1 + w ' ^ ^ \ + w' \ + w ^ ' 

the estimate of the curvature term in (23d) is not included in these conditions; we refer to Corollary 6.4. 

The analysis of (23e) requires a thorough investigation of the terms on the r.h.s. of (7a). In this process, 
we may replace, in each term, the variables (^ab, ff"'', tr fc, cr"f,, /n, Ua) by the solution of the asymptotic 
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system ( '^gab, *'g°'',tr %, ^c^t, ^^^a)- To see that this replacement is possible, consider, e.g., the term 
^u"'dtUa of (7a). Inserting (22) the term iJLu°-dtUa becomes 

( V + t/5-4/3(l+-)^) ( yb ^ ^a, -a^ O^bc) ^ 0^^ ^ ^2«,/(l+^)+a„ 0^^ ^ ^2»/(i+^)+„„ _i j^^^^^ ^ ^^^^j^ 

Because of (24e), the power / goes to zero as t faster than °Ua; the same is true for 
£2w /{i+w)+ai,-i comparison with dt and ^^-4/3(1+™) comparison with "z^. We conclude that the 
term ii+4/3(i+tu)-/} ^u°-dtUa of (23e) is of the required form J- +Q tdtU, where and Q converge to zero 
as 0{t^) for some positive power 5, if ^1+4/3(1+™)-/? o^o^a^^o^^ g q^^s^ 

We find that the term on the r.h.s. of (23e) with the slowest convergence is ^1+4/3(1+™)-^^^! _(_ ^u"". 
This term thus yields the most restrictive condition on which reads 

2 u) - 1 

/3 < 1 - = . 28 

Finally, in order for the r.h.s. of (23f) to be of the required form, we must impose the additional conditions 

2 „ 3u; + 5 ^ , , „ 10 - 1 

a a. < 1 + +fi= + /3 b a,„ < 3 - = ^ ^ • 29 

3(1 + w) 3(1 + w) 3(1 + w) 3(1 + w) 

Equations (23g) and (23h) yield a\> Q and a\ > 0. 

In the remainder of this subsection we consider the most complicated expressions of the r.h.s. of (23), 
which are the curvature terms in (23d). The estimates of the curvature terms are presented in a succession 
of lemmas. Throughout we use a coordinate frame to perform the computations. 

Lemma 6.1. 

°9ab, \b,c, \b,cd e O(tTO)), V', Y\c e O(rTO)). (30) 

Proof. The statements follow directly from the series expansions equations (19c), (19f) and the smoothness 
of the initial data. □ 

Lemma 6.2. 

°Rab e 0(1). (31) 

Proof. In a coordinate frame we have 

2"r%e = Y" ( V,c + V,6 - V.d) e o(i) , 



0(t 3(1+™)) =:0r„i,o6O(t3(l+"')) 
2°r°bc,d = \"^,d "Tebc + Y^ '^ebc.d G 0(1), 



0{t 3(i + ™))0(t3(i+™)) 0(t 3(i + ™))o(t3(i+™)) 

which gives = °r^a6.c - °r^,,,, + °r^d,or'i,b - or^db^r^^^, e o(i). □ 

Lemma 6.3. // a-y > ax > 0, and a-y > a\ > then 

^Rab e 0(l). (32) 

Proof. Assume a-. > a\ > 0, and > a\ > 0. We have 



s 



lab = Vb + V-( Vc7'b + Vl'a) e 0( "gab) C O(tW^), 
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Sgab ^ - V',c + V\c + fd Y\c) + it"--"^ (^rfc V + A'rfc V") e 0(i" W^), 

9ab,cd= gab,cd + -jt ^{ 9ae,cdl b+ 9be,cdl aj+I ^ { 5aeA be + ff&eA acj 

+ °5aeA%c,d + °9beWc,d) G O(i^) . 

The result then follows as in the proof of Lemma 6.2. □ 

Corollary 6.4. The curvature ierm i'^i? (5°^) on the r.h.s. of equation (23d) vanishes 

with t if OL^ is any positive constant and the conditions of Lemma 6.3 are met. 

It remains to collect conditions (24)-(29) and the result of Corollary 6.4: The reduced system (23) is a 
Fuchsian system of the form (21) that satisfies the requirements on the matrix and the functions if the 
constants (a-y, a-y, olu[= /?], au, cka, a\) are chosen according to the conditions 

4 _ 4 

< a\ < , < Q!A < a-, , + a,; < 2 — — r , + a,; < 2 — 7— , 

3(1 + w) ' 3(1 + t«) 

2 2 2 

/3 < 1 - — — , a, + p<l- — — , a„ < 1 + — — + p . 

l + W l + W 3(1 + ill) 

Under these conditions, there exists, for each choice of solution "VV = ( H7Q(,,tr "^fc, "cr^^, "yu, %a) of the 
asymptotic system, a unique solution U of the reduced system that vanishes as t — > 0. 

To obtain the statement (20) of Theorem 4.1 we set 

4 2 e 4 

a =a^ = 2-—— r-2£, = = /3 = 1 - — e, = - , a„ = 2-— — --2e(34) 

3(1 + w) l + W 2 3(1 + w) 

which is compatible with (33) for arbitrarily small e > 0, and insert these constants into (22). 
6.2 Symmetry of the solution of the Fuchsian system 

In the second part of the proof of Theorem 4.1 we need to come back to the problem of symmetry: 
Using the reduced system and its unique solutions that vanish as t — > 0, we have constructed solutions 
W = (^ab, tr fc, c";,, /Lt, ita) of the Einstein-Euler evolution equations (4) and (7) with modified curvature 
terms, and ^R, cf. the remarks following (22). A priori, it is not obvious that the constructed 

tensors gab and aab are symmetric tensors; however, if they are, then the solutions W we have obtained 
are solutions of the proper Einstein-Euler evolution equations. (The Einstein-Euler constraint equations 
are analyzed in the subsequent subsection.) 

Let °yV = ( "gab, tr "V^^, %„) be a solution of the asymptotic system as given by (14). Consider 
the unique solution U of the reduced system (23) that vanishes with t and define, through (22), the 
quantities W = (f^ab, tr fc, cr^jj, /i, Ua) . By construction, W satisfies the Einstein-Euler evolution equations 
with modified curvature terms. These equations imply equations for the antisymmetric parts of 9ab, g"'^, 
and aab. We have 

dtg[ab] = 2CT[Qb] - f (trfc)5[„b] , (35a) 
a^^I^^l = -2a[^5^^1 + |(tr k)g'°'^'\ , (35b) 

dt(y[ab\ =■ |(tr/s)CT[Q6] + 2ayac(y''b\ + |(w + l)/^'«^ff[ab] + 9{ab\- (35c) 



12 



The terms cr^'^cg'^''^ and a[ac<^'^b] in (35b) and (35c) can be expanded by splitting both aat and g"'' into 
its symmetric and antisymmetric part. Each term in this expansion then contains at least one of the 
anti-symmetric tensors crj^j^] or gl"''!. We have 

^[a^gCb] ^ gi-c)^^^^^g[db] ^ g[-c]^^^^^g(db) ^ g[ac]^^^^^gm + gi^^) g^ ^^^^^ ^ (35d) 
'^[ac'^'^b] = 0'(ac)5^''''^0-[d&] + '^[ac]9^'"^^ (^{db) + Cr[ac]9^'"^^ (^[db] + '^{ac)'^{bd}9^'"^^ ■ (35e) 

Let g, g, and a be constants and define the antisymmetric tensors fi^^ := tSg[ab], ^g'' '■= t^g^°''^\ and 
fJ'^jj i'^cr[a6]- From (35) we obtain a system of equations for these tensors, 

td,nl, + {-g - Ij^PI, = 2t^+'^-^ni, [l+'-^t trk\n\, , (36a) 

td.ni" + (-5 + = -2ti+V^/''l + ^ [1 + i±H t tr fc] , (36b) 

tdtK, + (-a + 1^)1]^, = 2ii+V[,£a^,] + 3^[1 + ^ttrk]n:, 

+ ^t'+'^-^ ni, - 1(1 + wWt'+'^-<'ni, , (36c) 

where we use (35d) and (35e) to express t-'^+s o-[a^^cb] aj^^j t^^'^o-[acO''^b] in terms of ri^^, ^g'', and 51"^^, 
according to 



By an appropriate choice of the constants g, g, and a we are able to achieve that the system (36) is of 
the Fuchsian form (21), i.e., (i) the factors in front of ri^^, ^g''' ^ab '^^ ^^'^ ^-^^ sides of (36) are 
positive, and (ii) the factors of ^g'', and fi'^j, on the r.h. side vanish like a power of t with t — > 0. 
(As we have established, the leading asymptotic behavior as i — > of the terms on the r.h.s. of (36) is 
given by replacing g(ab} by "ffab and using (19), and analogously for the remaining expressions.) From 
the l.h.s. of (36) we have 

4 4 2 

3(1 -fw) 3{l + w) 3(1 + w) 

From the r.h.s. of equations (36a)-(36e) we obtain the conditions 

l + g-a>0, l-a-5>0, 1 + g-a- ^ , >0, (37b) 

3(1 + wj 

^ + ^-9 - on 1 ^ > > 1 + a-g- ^ > . (37c) 
3(1 + iy) 3(1 + w) 

The third set of conditions we impose on the constants g, g, and a is that r2^f,, ^g'', and fi'^^ vanish in 
the limit t — > 0. E.g., from (19) and (20) we have 

^Ib = %cl\ e 0(t«+2-2e) ^ 

while it turns out that Q^'' e 0{t^+^^n^~'^'') and fl^f, e 0(1"+^-^/'^). This gives 

5 + 2-2e>0, g + 2- ^ - 2e > , a+l-|>0. (37d) 

3(1 + w) 2 
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The positive constant e can always be chosen smaU enough so that it is possible to satisfy the collection 
of conditions (37) for an arbitrary w > 1. An explicit choice is 



Sui _ 1 — w 1 + 3w 

a — 



" 3(1 + w) 1 + w' 6(1 + 

Since (36) is a Fuchsian system of equations for the variables fi^^, ^g' ■, ^"ab- ^^^^ solution fi^^ = t^g[ab]; 

a[ab] we have in our hands, is the unique solution of (36) that vanishes as t 0. 
However, the trivial solution is an explicit solution of (36) with the obvious property that it vanishes 
in the limit t — >■ 0. The uniqueness result thus ensures that Jl^j^ = t^g\ab\ = 0, fig'' = V^g^""^^ = 0, and 
n^fj = ^^[ab] = 0. We thus obtain symmetry of the metric and the extrinsic curvature, which are thus 
solution of the proper Einstein-Euler evolution equations. 

Remark. The proof of the symmetry of the constructed solution that is given here is inspired by the 
considerations of [3] for the stiff fluid case, i.e., w = 1. However, the relevant arguments in [3] are 
incorrect, because the equations in [3], which are a variant of (35), lack a term that is quadratic in the 
extrinsic curvature. This problem can apparently be fixed, see [4]; we note, however, that the arguments 
we follow here for w > 1 fail for w = 1. The conditions (37) lead to a contradiction for w = 1, because 
in that case we find 

/ 8 , 1 _ 4 1 

g-a> -{I - —— rj = q , g - a < I - —— r = -. 

3(1 + 10) 3 3(l + iy) 3 

In addition, the first term on the r.h.s. of (36d) and (36e) behaves like t("'~i)/(i+™); this function vanishes 
for t — )• in the case of an ultrastiff fluid, i.e., w > 1, but does not vanish when w = 1. 



6.3 The constraints 

To conclude the proof of Theorem 4.1 we consider the Einstein-Euler constraints (5). We need to show 
that the constructed solutions W = {(jab: tr fc, a^f,, /x, Ua) of the Einstein-Euler evolution equations satisfy 
the Einstein-Euler constraints (5). If this is the case, then the solutions W we have obtained are solutions 
of the Einstein-Euler equations (4)-(7) and Theorem 4.1 is proved. 

Consider W ~ {gab, tr k, cr°j, Ua) and let 

-a\a\ + j{trk)^ + R-2p, = _Vba''„ - |V„ tr fc - , (38) 

where p and ja are given by (6). The Einstein-Euler evolution equations imply equations for these 
constraint quantities. We have 

at$ = 2trA:$ + V"4'a, ft^a - tr fc *a - 5 V^* . (39) 

Let ip and ip be constants and define fi* ~ t'^<I> and fi* = t'^''^a- From (39) we obtain the system 

tdtn-" + - (^) 17* - ^ [1 + i±H t tr fc] f}* + t'+'^-^V-nl , (40a) 

td,nl + i^pl ^ ^ [1 + i±Ht trfcjf]* - , (40b) 

where V°f2* can be expressed as 

= {dag''')nY + g'^'dant + y'"'g''{^agbc)n^ ■ (41) 

It is straightforward to show that the system (40) is Fuchsian when the constants (p and ip satisfy a number 
of conditions. (These conditions follow from (40) in the usual manner; in (41) we use the constructed 
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behavior of the metric and its spatial derivatives, see (20).) In particular, there exist appropriate choices 
of and ip that are compatible with the additional conditions that fi* = and 51* = t'^'i/a vanish in 
the hmit t ^ 0. To obtain these conditions we note that -°cr"^ V'^ + |(tr °fc)2 - 2 "/i = 0, cf. (10), which 
results in 

n'^ =t^{- a\a\ + |(tr kf + R- 2p) 

where "higher order in t" denotes terms that converge to zero faster than at least one of the terms that 
are given explicitly. From (19) we see that it is the curvature term that determines the leading behavior 
of rj* in t. We obtain fl* e (9(^v-4/3(i+»))_ Analogously, we have 

^t^{-{^b- Vb) "a\ - tl-4/3(l + .«)-s/2 O^^^b^ _ 2^2-10/3(l + u,)-ey^^ 



- (1 + w) % (Vl + u2 - 1) V - (1 + u,)t4-10/3(l+«;)-2e 

- (1 + u;)t^-i°/3(^+"')--V V + higher order int). 

The leading order comes from the term V^^^^. We thus obtain 17* e 0(t'A+(3»"-i)/3(i+u>)-5/2^^ rpj^^ 
required conditions are satisfied if the constants ip and -0 obey the inequalities 

4 4 12 , 4 , 2 6 , , 



3(1 + w) l + w 3{l + w)' ^ 3(1 + w) 2^ ' l + w 3(1 + w)' 

l + ^-i'- TrTT^ 1 + ^ - > . (42b) 

3(1 + wj 



Assuming that e is small enough (which we can always do) we are able to ensure that the required 
conditions are satisfied; an explicit choice of constants is (p = ip = 5/3(1 + w). Then (40) is a Fuchsian 
system of equations for the variables £7*, il*, and O* = and il* = t^'ia represent a solution of 
that system that vanishes as t — )• 0. However, since the trivial solution is an explicit solution of (40), the 
uniqueness result guarantees that f2* = t"^^ = and 17* — t'^'i^a = 0, i.e., the constraints $ = and 
= of the Einstein-Euler system are identically satisfied. 

The solution W = [gab,iJ^ k, a°-f^, fi,Ua) of (20) is thus a solution of the Einstein-Euler equations (4)-(7) 
and Theorem 4.1 is proved. 
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